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Abstract 

In a recent work we have discussed how kinetic theory, the statistics of classical particles obeying 
Newtonian dynamics, can be formulated as a field theory. The field theory can be organized to 
produce a self-consistent perturbation theory expansion in an effective interaction potential. In 
the present work we use this development for investigating ergodic-nonergodic (ENE) transitions 
in dense fluids. The theory is developed in terms of a core problem spanned by the variables p, 
the number density, and B, a response density. We set up the perturbation theory expansion for 
studying the self-consistent model which gives rise to a ENE transition. Our main result is that 
the low-frequency dynamics near the ENE transition is the same for Smoluchowski and Newtonian 
dynamics. This is true despite the fact that term by term in a density expansion the results for 
the two dynamics are fundamentally different. 







I. INTRODUCTION 



In previous work jlj], referred to as ND1 from hereon, we showed how the Newtonian dy- 
namics (ND) of simple fluids can be treated using field-theoretic methods 2|. Here we apply 
these methods to systems in equilibrium with an emphasis on the ergo dic-nonergo die (ENE) 
transition which is generally associated with the mode coupling models {3 5] for glassy dy- 
namics. We analyze the long-time behavior of the density- density time correlation function. 
For high densities we find that there is a slow dynamical regime which is very similar to 
that found in the case of Smoulokowski dynamics (SD) . In the shorter-time regimes 
the SD dynamics is considerably simpler than ND. The complications in the ND case arise 
because of the unwieldy forms of the three-point vertices even for the ideal-gas limit. In 
the long-time limit we find that the three-point cumulants and three-point vertices satisfy 
fluctuation-dissipation relations which are very helpful in treating the late-stage kinetics at 
high densities. 

After determining the nature of the fluctuation dissipation relations j§| for the two-point 
cumulants and vertices, we derive an exact expression (kinetic equation) satisfied by the 



density- density correlation function. The form of this kinetic equation is the same as found 
in projection operator treatments^] of the problem. This approach involves the introduction 
of memory functions. The significant difference between our work here and the projection 
operator method is that we obtain expressions for the memory functions which are suitable 
for perturbation theory analysis. It is not expressed in term of "projected dynamics". 
Instead the memory function is expressed as the "BB" component of a two-point vertex 
function Tbb associated with what we call a core problem. The core problem involves the 
kinetics of the particle density p and a response field " B" . 

We then carry out perturbation theory in terms of a pseudo potential and find at second 
order that the memory function kernal Kbb in the low-frequency regime is proportional to 
a quadratic-form in the full density- density correlation function. This suggests a non-linear 
feedback loop similar to that found in mode-coupling theory 10|, QjJ and in treating SD 6] . 

When we focus on long-time solutions we find that the developments for SD and ND are 
the same. In SDENE 6] we analyzed the leading contribution near an ENE transition. In 



SM 



12] it is shown the ENE-transition is accompanied with a two-step decay process. SM 



shows that there are a pair of exponents, a and b, which characterize the slow-time power- 



1 



law solutions. They treat in detail the case of hard spheres at second order in perturbation 
theory. We find here that SD and ND systems share the same ENE transition according to 
our results at second order in perturbation theory. This point is amplified in Ref,. 



13 



where 



higher order terms in perturbation theory are treated. We speculate that there is a whole 
class of systems which share the same statics and ENE transition. 

II. FIELD-THEORETICAL FORMULATION 
A. Newtonian Equations of Motion 

We discuss the kinetics satisfied by a system of N particles with mass m, position Ri and 
momentum Pi we consider Newtonian dynamics, with equations of motion 

mRi = P l (2.1) 

Pi = fi , (2-2) 

fi = -ViU , (2.3) 
U = ^V(R i -R j ) (2.4) 

i,3 



where the force is given by 



the total potential energy is 



where V(Ri — Rj) is a standard pair potential between the z-th and j-th particles. Starting 
with these equations of motion we can reexpress the problem in terms of a path-integral 
formulation. 

B. Path- Integral Formulation 

In FTSPD and ND1 we introduced a field-theoretical formulation for classical many- 
particles dynamics. Here we summarize the results for Newtonian dynamics. The grand 
partition function for a core set of dynamical fields = $ is given by 

Z ? = E Z m TTe ~ A+H '* ( 2 - 5 ) 

N=0 
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where the trace is over the set of phase-space variables Ri(t) and Pi(t), and the conjugate 
MSr[i^] variables Ri(t) and Pi(t). z is the fugacity and Hi is an external time-dependent 
field that couples to the set of dynamical variables $ a . The MSR action A is given by 

i J 

In the present work the set {3> a } stands for the core variables {p, B} which are collective 
properties and are expressed as a sum of one point functions. Explicitly the particle density 
is given by 

TV 

p(l) = ^5(a; 1 - J R J (t 1 )), (2.7) 
i=i 

and the response field B: 

B(x,t)=i^P i .-—6(x-R i (t)) . (2.8) 

The -B-field is important in the present analysis and is somewhat unfamiliar. We note that 
B is the longitudinal component of the vector field g(x,t). The latter is expressed as a 
collective density corresponding to the hatted MSR field Pj (with a factor i ) conjugate to 
the momentum Pj of the i-ih particle. 

N 

g(x,t) = Y,^(x-Ri(t)) (2.9) 

i=i 

This g(x, t) however is not the field conjugate to momentum density g(x, t) used in the usual 
MSR formulation of equations of nonlinear fluctuating hydrodynamics. The hatted fields in 
the MSR field theory are generally viewed as convenient mathematical tools for ensuring the 
that the stochastic equations of the dynamics are satisfied in the field theory. These auxiliary 
fields give rise in a natural way to a set response functions which are related to correlation 
functions through suitable fluctuation dissipation relations. From such considerations, the B 
is best interpreted here as a suitable mean to obtain linear fluctuation dissipation relations 
facilitating the analysis of the dynamics developed here. 

In general it is assumed that we can express $ Q as a sum of one particle contributions 



iRi{t)\Ri{t) 



P(t) 



iP t {t)\Ht)- m 



dt 



(2.6) 
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iV 



i=l 

where with the roman index i we indicate that the corresponding quantity is a single-particle 
property. In the above equation we have used a compact notation where the index a labels 
space, time and fields p or B. We maintain the notation from here on. For a = p we obtain 
the collective density for which we have 

<j>®(l) = 8(x 1 -R i (t 1 )) . (2.11) 



For a = B field. 
For the SD case [2 



p 



;he corresponding </>g (1) is strongly dependent on the type of dynamics. 



while for the ND 



= DiiR^V^ + ~V|^(xi - (2.12) 
l| case the corresponding result is 



0f (1) = -iA(*i)V ccl 5(a; 1 - Riih)) (2.13) 

where -Pi(ti) is the MSR conjugate momentum for the selected particle. 

For our purposes the MSR action A defined in eqn. (I2.6P can be written as a sum of two 
parts, 

A = Ao + Ai . (2.14) 

The first term on the RHS Ao is the noninteracting MSR action. The interacting part of 
the action is given by 

Ai = \Y<J dld2$ Q (lK,(12)<M2) (2.15) 

where the Greek labels range over p and B. The interaction matrix, for systems in equilib- 
rium in the distant past, is defined just as in SD, by 

<7 Q/3 (12) = V(12) [5 aB 5 Pp + S^Spa] (2.16) 

and 
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V{12) = V(xi-x 2 )S(t 1 -t 2 ) ■ (2.17) 

Notice that the response field B is chosen such that the interaction part of the action has 
the form given by Eq. (12.15j) . 

For studying the dynamics of a many-particle system, we will be interested in the cumulants 
generated by corresponding generating functional 

W[H] = \nZ T [H] . (2.18) 

The one-point average of a field $ Q defined in eqn. (12.101) in an external field is obtained in 
terms of the functional derivative 

G a = i^ a ) = ^ f W[H] . (2.19) 
The full cumulants G a p, G^ etc. are defined by successive functional derivatives of G a , 



G a p 
G a f3~f 
G a 0...8 



5 
5 

5IL 



G, 



G 



a/3 



7 

_5 5_ 

W a Wp 



We also need to deal with the single-particle 



5H, 



■W[H] 



(2.20) 
(2.21) 
(2.22) 



15| quantities 



N 



J a rp ■ ■ ■ rs 



(2.23) 



i=l 



III. FLUCTUATION DISSIPATION RELATIONS 



In this section we explore the nature of fluctuation-dissipation relations in case of New- 
tonian dynamics assumed to be in equilibrium. The correlation functions between the fields 
in the filed theory is primarily determined with the corresponding action functional of the 
associated field theory. The MSR action without initial conditions and external fields is 
given by Eq. (I2.6j) . We now explore the invariance of the action A under certain symmetry 
operations. In a stationary state we have time-translational invariance and time-reversal 
symmetry. 



A. Symmetry transformations 



1. The time-reversal symmetry T is defined as, 



TRi{t) 
TRi(t) 
TP t (t) 
TPi{t) 



Ri(-t) 

-Ri(-t) 
-Pi(-t) 

H-t) . 



(3.1) 



We consider how the MSR-action A changes under time-reversal 



x = ta = J2 j dt i {-M-t)} { 



dRj(-t) | Pjj-t) 
d(—t) m 



(3.2) 



Letting t — > — t in the time integral, we obtain A' = TA = A if the limits of integration are 
symmetric. This includes t2 — > oo and t\ — > — oo here. We conclude that the MSR action 
remains invariant under time-reversal. 

2. In ND1 we have discussed the Fluctuation-Dissipation symmetry (FDS) for the field 
;heory corresponding to the MSR action functional defined in Eq. f)2.6p . It was demonstrated 
1(J that the following transformation r keeps the MSR action invariant : 



rRi(t) = Ri(-t) (3.3) 

rRi(t) = -Ri(-t)-0fi(-t) (3.4) 

rPi{t) = -Pi{-t) (3.5) 

rPi(t) = P t (-t)-t^P t (-t) . (3.6) 
m 



We have therefore tA = A. 

3. The inversion symmetry denoted by T\ is defined as 
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nRiit) = -Riit) 
nPi{t) = -Pi(t) 

nP(t) = -P t (t) . (3.7) 

The following invariant properties of the Action A easily follows : T\A = A. 
4. Finally we define the symmetry r 2 as 



T 2 Ri(t) = Ri(t) 
r 2 Ri{t) = -Ri(t) 

T 2 Pi(t) = P(t) 

r 2 P(t) = -Pi(t) . (3.8) 
The following invariant properties of the Action t 2 A* = A is satisfied. 

B. Fluctuation-Dissipation symmetry 

From the definitions ( 12.7ft and ( 12. 8ft it follows that under the transformation r the core 
fields p(x,t) and B(x,t) change as 



rp{x,t) = p(x,-t) (3.9) 
tB(x, t) = B(x, -t) - (3.10) 

In ND1 we had demonstrated the following Fluctuation dissipation theorem involving two- 
point correlation functions between the p and B fields. 

G fB (t - = 9(t - t')(3^G fp (t - t'), (3.11) 

for any function f[p\. 

Here we are interested in higher-order FDR which are most easily expressed in terms of 
Fourier transforms. Consider the mixed correlation function 
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C BB ... Bp ... p (12...££ + l...n) = (B{l)B{2)...B{l)p{t + l)...p{n)) 

where 



(3.12) 



N 

p(j) = p{q j ,u j )= / dte^Y^e 
J i=i 



-iqj-Ri(t) 



and 



i=l 



-iqj-Ri{t) 



Under the transformation r we have 



(3.13) 



(3.14) 



tB(\) = B(qx, +i^uxp(q 1 , -Ux) 
Hence the mixed correlation defined above in Eq. (I3.12p is obtained as 



(3.15) 
(3.16) 



C BB ... Bp ... p (12...££ + l...n) = {tB(1)tB(2)...tB(£)tp(£ + l)...rp{n) 

'[B(l)+iPu lP (i))[B(2)+il3u 2 p(2)) 
\B{£) + tPu eP (£)]p(£ + l)...p(n)) (3.17) 



where j = qj,—Uj. Multiplying these out, each correlation funcion has arguements with 
tildes. Using the result 

nr 2 A* = A (3.18) 

and 

r 1 r 2 B*(l) = B(l) (3.19) 

7iT2P*(l)=/0(l) (3.20) 

we find the simple result 



C BB ... Bp ... p (12...££+l...n) = ([B(l)-iPu 1 p(l))[B(2)-0u 2 p(2)) 



.. [B{£)-if3uj e p{£)}...p{£ + l)...p{n) 



(3.21) 



If we define 



= B(i) - iPu iP (i) , (3.22) 
the mixed correlation function is obtained as 

C BB ... Bpp {\2 ...£,£ + 1 ... n) = (^(1)^(2) . . . tf; (£)p(£ + 1) . . . p{n)) (3.23) 
One easily sees that this holds for the cumulantsQ] 

G BB ... Bpp (12 ...£,£ + 1 ... n) = G^„.„ p (12 ...£,£ + 1 .. .n) (3.24) 
and for the single-particle quantities 

Gbb...b pp (12 ...£,£+ 1 ... n) = ^...^... p (12 ...£,£ + 1 ... n) . (3.25) 

Two-point cumulants 
Using the FD relations at the two-point level we find 

G Bp {12) = G* Bp {l2) + ifaGTJU) . (3.26) 
Due to translational invariance in space and time 

G af) (12) = G a p(l)5(l + 2) (3.27) 
G pB (l) = G* Bp (l) . (3.28) 

We see that Eq. (l3.26p can be written as 

G Bp (l) = G pB (l)+0cJiG* p (l) . (3.29) 
Since G PP {1) is real, using eqns. (13. 28ft and (13. 29ft . we have the conventional result 

ImG Bp (l) = -ImG pB (l) = ^G pp (l) (3.30) 
which can be used to construct the dispersion relation 



p v ; ' 2tt ui - Co + ir] y ' 



Taking the inverse time Fourier transform we find 

GpB(q,t)=e(t)P^G pp (q,t) . (3.32) 
Using the FD relation (13 .31 p we obtain for G p s(g, 0) 

G pB (q,w = 0) = / ^— = -PpS{q) 3.33 

J Z7T —UJ + IT} 

where S(q) is the static structure factor at wave vector q and p is the average density so that 
the RHS represents the correlation of density fluctuations at equal time. From eqn. (I3.29p . 
we obtain the very useful result 

i{3cuG pp (q, u) + G pB (q, u) - G Bp (q, u) = (3.34) 

Eqn. (I3.34p is written in a compact form as J-[G] = 0, where the functional T acting on the 
elements of a 2 x 2 matrix A pu is defined as follows 

J r [A] = Ti[(i/3ujl + e)-A(q,uj)} . (3.35) 
We have introduced the traceless anti-symmetric matrix e in the RHS of eqn. (fdr-func) : 

Bb p = 1, £ p b = -1, £ MAt = 0. (3.36) 
X is the 2x2 identity matrix. The functional T is defined as 

f[A\ = Tr [(-z/3wX + e) -A(q,u)] (3.37) 

where the transpose matrix e pu = e up . Vanishing of the above defined functional J 7 [A] then 
ensures the fluctuation-dissipation relation (FDR) among the of the elements of the matrix 
A. It is straightforward to show that the relation 

F{A] = (3.38) 

is satisfied for the two-point propagator A belonging to any member of the set {G^°\ G,Q}. 
With the condition ( I3.38|) the elements of A will be referred to as FDR matrix propagator 
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(FDRMP). We show in App endix ICl how to build new propagators which also obey the same 
FDR. If A a p(q, u) and D a p(q,u) are respectively two FDRMP then the matrix H a p defined 
in terms of the interaction matrix er,,„ as 



is also a FDRMP since it satisfies the relation JF[H a p\ = . This property of composite 
propagators following the same Fluctuation-dissipation relation is an important ingredient 
in the analysis of memory functions discussed later in this work. 

2. Two-point Vertices 

The two-point vertices T a p(l) are defined as the corresponding elements of the inverse of 
the general two-point correlation matrix G a p : 



For a = p and (3 = B, using Gbb(X) — 0[17|, it follows from the above equation that 
r pp (l) = 0. For other choices of a, (3 we obtain, 




(3.39) 




(3.40) 



T pB {l)G Bp {l) 
T Bp {l)G pB {\) 
T BB (l)G Bp (l)+T Bp (l)G pp (l) 







1 



1 



(3.41) 



(3.42) 



(3.43) 



Using Eqs. (I3.4ip and (I3.42p . Eq. ( I3.43|) can be written in the form 



G PP (1) 



G pB (l)T BB (l)G Bp (l) 



(3.44) 



We easily find 



T Bp (q,u) = T* pB (q,u) 



(3.45) 



Starting with Eq. (l3~26|) : 
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G Bp {l) = G Bp {l)+iPuj l G; p {l) (3.46) 
we can use Eq. (l3T4ll . (l3~42l) and ( l3~45l) to write 

' ' +i/3 Wl G w (l) (3.47) 



r>(i) v Bp {\) 

Furthermore using Eq. ( 13.43]) and Eq. fl3.47p . and canceling a common denominator, gives 

iPuj l T BB {l) = T pB {l)-T Bp {l) . (3.48) 
In terms of the functional T defined in Eq. f!3.35j) the above relation is obtained as 

•F[T] = , (3.49) 
where the matrix vertex T a p is defined in the {p, B} domain. 

3. FD symmetry and Three-point cumulants 

There are substantial differences in specific cumulants for the different realizations of 
simple-fluid dynamics. For example the noninteracting three-point cumulant G^ p (123) for 
SD and ND are quite different. Despite this fact we find that both sets of cumulants satisfy 
the same fluctuation dissipation relations. 

In treating the three-point cumulants involving the fields at three different points, we 
adopt the following notation (in Fourier space) 



G PrPoPk = G p PP (ijk) = G ppp (jki) = G ppp (kji) (3.50) 
G PiPjBk = G ppB (ijk) = G pBp (jki) = G Bpp (kij), (3.51) 
G BlBjBk = G BBB (ijk) = G BBB (jki) etc. (3.52) 

where the set {ijk} =>- {123}, {231}, and {321}. The factor (3 will for now be absorbed in 
the frequency u. For cumulants involving three B fields we obtain using the basic fluctuation- 
dissipation relation (13.261) the following relations : 
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G Bl B 3 B k - G* B . BjBh + iu)iG* p . B . Bk + iujG* B . pjBk + iu k G* BiBjPk (3.53) 
- ujiUjG*.^ - u) j u k G* B . p . pk - u k UiG* kB . B . - iuiUjjUkGp.p.^ 

Since G BBB = G* BBB = 0, using this in the complex conjugate of Eqn. (I3.53P gives the 
relation 

= iuJiG PiBjBk + iUjG BiP:jBk + iu k G BiBjPk + ^,^/ : l . l ,,n l (3.54) 
+ UjU k G BiPjPk + bJ k bJiG PkBjPi - ibj i bJ i bj k G PiPiPk 

Similarly for the cumulants respectively involving two and one B fields are obtained as 
follows : 



G Bl B ]Pk - G* B . B . ph + iWiG*. B . pk + iujG* B . p . pk - Uiufi*^.^ (3.55) 

G BiPjPk = G* Bip . ph + iuiG* pip . pk . (3.56) 

To simplify the above relations involving the three-point cumulants, we define /ij = 
iUiG PtBjBk , Oi = UjUJ k G BiP] p k , and v = UiUjUkG PlPjPk . In terms of the newly defined quanti- 
ties we rewrite the Eqs. f l3.55p - fl3.54p respectively as 

— nl = fi k + (Tj + ai-iv (3.57) 

ai = o* + iv (3.58) 

/ij + fij + + (?i + Uj + a k = iv (3.59) 

where v is real. Solving the above equations we obtain 



a i = u/2 (3.60) 
+ 
~2 



»i = — —z — ( 3 - 61 ) 



& = \ (3.62) 
where \it = \i% + iij + fi k is the sum of the /i's. These results seem rather general. 
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4- FD symmetry and Three-point Cumulants 

Next we analyze the relations between the different three-point vertex functions. For this 
we use the general definition of the three-point cumulants and the corresponding three-point 
vertex function as, 

G aiHVk = —G aa i {i)G^> {j)G uu > {k)Y a i ilx i. v i k (3.63) 

where «j stands for the field a at the point i and so on. The fields a,p and v are respectively 
taken from the set {p, B}. In the RHS of Eq. ( I3.63|) we have used the notation G pu {i) 
following Eq. (I3.26p . i.e., G py {ij) = G pv {i)5{i + j). Since Gbb = 0, taking all the three 
fields to be B at the respective points i, j, and k we obtain, , 

G Bt B 3 B k = -G Bp {i)G Bp (])G Bp (k)Y PiP]Pk (3.64) 

Since G BiBjBk = it follows from the last equation that T p . PjPk = 0. The expression ( I3.63P 
for the various three-point cumulants in terms of two-point cumulants simplifies due to the 
vanishing of T ppp . Let us first consider the expression (I3.54p involving three B fields. In 
the appendix [A] we show that by a) using the expression (I3.63j) for G apu in Eq. (I3.54p . and 
b) by collecting the various coefficients of the different vertices T's, we obtain the following 
relation between the three-point vertices. 

i(uiT BtPjPk + UjT PiBjPk + LO k T PiPjB kj - iu i u j u k T BiBjBk 

+ {uiU j r BiB . pk +U j U k T pi B j B k +UkUiTB iPj B k } = Q ■ (3.65) 

Next, we consider the relations (I3.55P which contain three-point cumulants involving corre- 
lation of two B fields. 

G BlB]Pk = G B . iBjPk - iuiG PtBjPk - iujjG BiPjPk - UiUjG PiPjPk (3.66) 

Using the the three-point cumulants corresponding to this case in terms of vertex functions 
as defined in Eq. (I3.63P and collecting the coefficients of the various vertex functions, we 
obtain in Appendix |A] the result : 

iuj kT* PiPjBk + iUiT BiPjPk + iujjY Pi B jPk + u i u j T Bi B jPk = . (3.67) 
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By doing circular permutations of the indices {i, j, k} in Eq. (I3.67P we obtain the following 
relations involving the vertices with one I?-field leg: 



iuj i T kp jPk + iuj i T PiB jPh + iu k T p . p . Bh + uj j uj k Y PtBjBk = (3.68) 
iuj 3 V * Pl B jPk + i u k T PiPjBk + mT BiPjPk + u} k UiT BiPjBk = (3.69) 

The above relations are further simplified with the following notations : 



4>i = iuiT BiPjPk , (3.70) 

tpi = UjUj k T p%BjBk , (3.71) 

P = -iui i u j uj k T BiBjBk . (3.72) 

We obtain from Eq. (I3.65f) using the above notation the result 

P + 0t + ^t = O (3.73) 

where (fir = 4>i + 02 + 03 and ipx = ^i + "02 + "03- In the simplified notation the Eqs. (I3.67P 
- (I3.69P respectively reduces to the cyclic forms 



- 4>* + h + 4> k + = . (3.74) 

Looking at the imaginary parts of these equations we have ip" = — 0^, while for the real 
parts —§\ + 4>'j + <P'k + i'i = 0- The last equation has the simple solution 

4>' 3 = -\{^ + ^ k ) (3.75) 

In terms of the original notation, using the defining relation (I3.63|) for the different vertices, 
we obtain the following set of nonperturbative relations among the various vertex functions. 



UjUkF piB:jBh 



u%u k Y" BiPiBk 



UiUjF BiBjPk 



-(uiF BiP . Pk + Ujr' p . B . Pk + u) k r' PiP , Bh ) 



(3.76) 
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The vertex function T BiBjBk is real, i.e., T BiB , Bk = 0. The real part satisfies the relation 

^^kV' BiBjBh = -2[u i I J BiPsPh +u j r , „ B .^+u k I J K Pj B k ] ■ ( 3 - 77 ) 

IV. KINETIC EQUATIONS 

We now present starting from the basic equation (13.401) . the derivation of a kinetic equa- 
tion for time evolution of the density autocorrelation function G pp (q, t). This kinetic equation 
involves the different contributions to the vertex functions IY,-. The analysis is similar to the 

n n 3 

discussions in SDENE[6[ and MMS[7j. These developments start with a discussion of the 
form of T a p(q, u) in the time domain. If we look at 

r ^> - to • (41) 

The response function G p b vanish algebraically as u — > oo, and hence the two-point vertex 
Y Bp diverges in this limit. We assume that 

lim T Bp (q, u) = -D q co 2 + iuA q + C q (4.2) 
where the coefficients A, C and D are real and not universal. We define the "local" quantity 

r^(g, u) = -D q u 2 + iuA q + C q (4.3) 
for all frequencies. The corresponding subtracted quantities are obtained as 

T%{q ) uj)=T Bp {q ) uj)-T%{q ) uj) (4.4) 
which vanish for large frequencies. At low frequencies 



rSJfa u = 0) = T Bp (q, u = 0)- r®(q, u = 0) . (4.5) 

Using the Eqs. 1ET351) and (jOl . we obtain T Bp {q,0) = ~[pf3S(q)y\ From Eq. (@~3]) it 
follows that T Bp (q, 0) = C q and hence using Eq. (I4.4p we obtain 

T^ p (q,u J = 0) = -C q -[pPS(q)}- 1 . (4.6) 
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Next, assuming that the FDR holds locally, we obtain from Eq. (14. 3 j) and FD relation ( I3.48j) 



uT%(q,co) = -2coA q (4.7) 

Hence T^ B (q,0) = — 2/3~ 1 A 9 . Assuming A q > and taking inverse Fourier transforms we 
obtain in the time space the following results. 

rS,(g,f-0 = -p-'AAt-?) (4.8) 



d 2 d 
Dq W~ q dt + Cq 



5{t - 1') (4.9) 



where the T's in the time space are denoted with a hat. For the full two-point vertex 
functions Yb p and Ybb the following relations involving the local and subtracted parts, 

t Bp (q,t-t') = T%{q,t-t') + T^ p { q ,t-t') (4.10) 
t BB (q, t-t 1 ) = -r'A^t - + f W (g, t - t') . (4.11) 

It then follows for the subtracted parts that we have the dispersion relation ( since the 
FDR holds locally for Y%\ ) 

reflecting the fact that the vertex function Tb p is analytic in the upper half plane. From the 
above FDR it follows directly that 



rg(g,w = 0)=)8f«( gi 0) . (4.13) 
In the time domain the FDR (14.12ft reads 

fg(g,*-0 = W-f)^fW( g ,t-f) . (4.14) 

A. Memory Function Equation 

The density auto correlation function G pp is of particular interest in the present theoretical 
model for studying the slow dynamics of a dense liquid. The natural order parameter for 
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the ENE transition which is the primary focus of this paper is the long-time limit of this 
function. In this section we obtain an equation of motion for the density correlation function. 
First we express Eq. (I3.43|) in the time-domain 



rf„(<7, 0)G pp (g, t-t') + F% B (q, 0)G Bp (q, t-t') = * q (t, t') 



■<£) 



(4.15) 



where 



*(/. /') = - / dstW(t - s)G pp (s - t') - f dst%l{t - s)G pB (t' - s) . (4.16) 



In writing Eq. (I4.15P we have used the fact that the response functions are time ordered so 
that f Bp (t — s) ~ 9(t — s) and G p B(f — s) ~ 9(t' — s). We then use the fluctuation dissipation 
relations ( EHZ} and ( 14341) to obtain 



ds 



d_ 

ds 



0$ B {t - s) 



Gpp(s-t) 



dsf^ B (t-s)—f3G pp (t'-s) . (4.17) 

-oo C>S 



Integrating by parts [18( in the first integral gives 



¥(t,f) = f3T^ B (0)G pp (t-t' 



1 ds^ B (t-s)-^G pp (s-t') 



t' Q 

dspt { * B (t-s)—G pp (t'-s) 

-oo C ' S 



1 ds^ B {t-s)^-G pp {t'-s) 



where we have assumed that t > t' . We then have the kinetic equation 



rg(g,0) + /3r^(g,0) G pp (q,t-t')- I ds(3T%(q,t - s)^G pp (q,t' - s) = . 



=.(«) 



(4.18) 



ds 



(4.19) 



Since t > t' we have dropped the T B > B (q,0)GB P {q,t — t') term in the LHS of Eq. (14.151) . 
Using the relations f !4.13|) and (14. 6[) we obtain for the quantity within square brackets on 
the LHS of fl4TT9|) as 



rgte, 0) + 0*l(q, 0) = t%(q, 0) + rg(g, U = 0) . (4.20) 
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Using the relations ( 14. 6 j) and (14.91) respectively for the second and the first terms on the 
RHS, the kinetic equation then takes the form 



D 92 A 9 1 
9 dt 2 q dt pPS{q) 



G pp (q, t) + J dspt$ B (q, t - s)^G pp (q, t' - s) 







(4.21) 



We see that our dynamical problem is now in the form of a memory function equation and 

* is) 

the dynamic part of the memory function is given by the subtracted part T BB of the vertex 
function. The static structure factor is the same for all fluids with the same pair potential. 
The coefficients A q and D q can be constructed using perturbation theory. In the case of ND 
we have the ideal gas result: 



rg(..«> = -^ 



1 



where z = u/(V2qvo), Vq = y/kriT /m is the thermal velocity and 



(4.22) 



rx 

S(x) = 1 - 2xe~ x2 / due u2 - iy/ir, 
Jo 

In the large frequency limit 



nxe 



(4.23) 



i 



-f-) 2+i 



(4.24) 



and on comparing this relation with Eq. (14. 3 p we can identify at the lowest order, 



Mq) 

C (q) 



1 

Tp 



q 2 D (q) = /3mC (q) 



m 
P 



(4.25) 
(4.26) 

(4.27) 



where the equilibrium density is given by p. With the above identification we obtain the 
kinetic Eq. (I4.2ip in the form 



<9 2 

— +fl 2 



where we have defined 



d 

G pp (q,t)+ j dsT mc (q,t - s)—G pp (q,t' - s 







(4.28) 
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q (3mS(q) 

T mc (q,t) = -%^ B (q,t). (4.30) 



Note that this resembles the second-order oscillator equation of standard MCT [10|, III] 
without the bare dissipation coefficient. In the case of Newtonian dynamics the microscopic 
dynamics is reversible and the kinetic equation ( I4.28P has been obtained without projection 
in to a space of slow modes. Indeed as we discuss below in order for the ENE transition to 
take place we need to break the time-reversal symmetry. 

The kinetic equation is diagonalized using a Laplace transform defined as 

POO 

G pp (q,z) = -i dte lzt G pp (q,t) . (4.31) 
Jo 

The Laplace transformed kinetic equation (14.211) reduces to the form for the normalized 
density correlation function F(q, t) 

F( n 7 \ G pp {q, z) z + iq 2 r mc (q,z) 

^ ) " pS(q) ~ z[z+iq 2 T mc (q,z)]-Q 2 q 1 " ) 

where pS(q) in the denominator of the LHS above is the equal-time correlation function and 

normalize F(q,t) to unity. We obtain the following integral equation for F(q,z) 



ril '- l) " r> /0 ° dte^/3V mc (q,t) 



l-zF(q,z) W q (q)J 

= ipS(q) / dte izt p 2 T^ B (q,t) 
Jo 

= -pS{q)p 2 Y^ B {q,z) . (4.33) 

We now turn to the mechanism which produces large F(q, z) and T^(g, z) as z — > 0. This 
involves determining the memory function as a functional of F(q,t). We define F(q) as the 
long time limit of the normalized density correlation function F(q, t) to denote the so called 
non ergodicity parameter. 



lim z ^o{zF(q,z)} ->■ F(q) 
From Eq. f |4.33|) we obtain the following integral equation for F(q) 
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(4.34) 



1 - F(q) 
T BB (q) is a functional of F(q) 



FUi) ~ pS(q)P 2 r$ B (q) . (4.35) 



hm z ^ {-zT^ B (q, z)} T^ B (q) (4.36) 

and is obtained from the explicit evaluation of the subtracted part of the vertex function 
T BB in a perturbation series. Solution of the resulting self-consistent equation, Eq. fl4.35p . for 
F(q) corresponds to what we call the ENE problem. This now primarily involves expressing 
T BB {q, t) self consistently in terms of the correlation functions. The self-consistent treatment 
of the dynamics constitutes the basic feedback mechanism and gives rise to the slow dynamics 



characteristics of the mode coupling theories 



10 
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V. SELF-CONSISTENT DEVELOPMENT 



A. General 



In Ref. 1 we obtained that the canonical partition function for the N particle classical 
system following Newtonian dynamics (ND) in the form 

r N 

Z N [H,h,h} = ^V{^ i )V{i! i )V{^f ) )e' Ao - Al+H - 4 ' +h ^ +h ^ (5.1) 
^ i=i 

where {^>i, ^j} are respectively the microscopic variables for the 2-th particle (i = 1, N) 
and its hatted counterpart in a Martin Siggia Rose (MSR) field theoretic formulation of the 
problem. In the above expression for Zn, the superscript in in the RHS refers to 
the initial state of the system. The MSR action is written as a sum of two parts as stated 
in eqn. (I2.14p . Aq is the quadratic part of the action including the quadratic contribution 
to the initial probability distribution. Aj is the interaction part of the MSR action given 
in eqn. (12. 15j) . The set {hi, hi} denote the currents conjugate to ^} for i = 1,...,N 
respectively. The collective field $ Q is defined in eqn. 1 12. 10p and H a is current conjugate to 
The canonical partition function (15.1 ft can be written in the convenient form 



Z N = Tl W e -Ai+H-S 
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(5.2) 



where we have introduced the N particle average as, 



r N 



(5.3) 



The partition for the interacting system is obtained using the identity 
in terms of operator 

&r = \ fdl ! d2 ^2H a {l)a aP (12)H p {2) (5.5) 

where H a (l) = [5 / 5H a {l)}. We rewrite the partition function in a form that allows us to 
formally carry out the average in eqn. (15. 2p using the functional identity (15. jj) . 



The sum over the degrees of freedom in eqn. (12. lip factorizes into a product of sums over 
the degrees of freedom of each particle. Together, these observations lead to the following 
result for the noninteracting Partition function, 

4 0) = TrW e ™ = ( Zl f . (5.7) 
where Z\ is the single particle partition function. 

Z x = Tr (1) e^-* (5.8) 

Working in the grand canonical ensemble, the grand partition function for the interacting 
problem given by (12. 5p is obtained in the form 

oo N 

Z T [H, h,h] = J2 % Z ^ H ' h ' M = eAreW ° • ( 5 - 9 ) 

N=0 

in terms of the single particle function Wq. In eqn. (I5.9p . z = Pq is the fugacity. The 
cumulants of the fields $ a are generated by taking functional derivatives of the generating 
functional VK[.ff] = In Z? with respect the fields H introduced in eqn. (12.191) above. 
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It was shown in FTSPD that the one-point average G a in a field $ a satisfies the funda- 
mental identity 

G a = Ti<p a e H -' t)+AW W. (5.10) 

where Tr is the sum over all the degrees of freedom including the conjugate MSR degrees of 
freedom for a single particle. The functional AW is defined by 

AW[H]=W[H + F]-W[H] (5.11) 

where 

F a = ^<T a p<j>p . (5.12) 

P 

The dependence of the theory on the interaction potential is controlled by AW [if]. This 
quantity is expressed in a functional Taylor-series expansion in powers of the potential V as, 

= E F °sk w m + E jw^fl + • • ■ ^> 

a a/3 H 

with F a given by Eq.f l5.12]) . For systems fluctuating in equilibrium the force matrix is given 
in terms of fourier transforms by eqn. (I2.16P above, 

<7ap{q) = V(q)[5 ap 5p B + 8 a B8p p ] . (5.14) 

The result ( 15 . 10[) was established in FTSPD (Ref 2.) using functional methods and ND1 
(Ref. 1) using algebraic methods. 

We conveniently introduce the set of cumulants: 

G -^ = m.wi-m w[m (5 - 15) 

to obtain the functional Taylor series expansion (I5.13P as, 

AW[H] = £ F a G a + \F a FpG a p + ^FaFffF-yGvfr + ... (5.16) 

a a/3 afi'y 

Clearly, in this form AW is treated as a functional correlation functions G a , G a p, G a/ 3 7 , ... 
and irreducible vertex functions. 
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We established in FTSPD a dynamic generalization of the static Ornstein-Zernike 
relation 19] given by eqn.(59) there. Inserting Eq. (15.16j) in eqn. (15.151) we obtain 



Otla — 

' 7 

where Q a p is given by Eq. ( I2.23|) and is related to the self Is- correlation function: 

<^ = Tr^e^+ AW . (5.18) 



20], or dynamic direct correlation function in eqn. (I5.17P is obtained 



The memory function 

as 

c a/3 = Tr<p a e H ^ +AW 4rAW. (5.19) 

Since AW can be treated as a functional of Gi we see at this stage that we have available a 
self-consistent theory. If we define the matrix-inverses 



TgsGgp = 5 aP (5.20) 

.5 

^lasGsp = 5 a p (5.21) 

<5 

then the two-point vertex is given without approximation by 

T a p = 7a/3 + K aP (5.22) 

where 

K a p = - ^'jascgp (5.23) 

8 

is the collective contribution to the vertex function V a p. From eqn. (I5.19P it follows that we 
can obtain the c a p and hence K a p as an expansion in the interaction potential by expressing 
AW in a similar series 

AW = AW {1) + AW {2) + AW {3) + ■■■ (5.24) 
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B. Collective part of vertex function at First-order 

Going to first order in pertubation theory we have 

AW (1) = • ( 5 - 25 ) 

a 

Then we need to compute 

fL 1] = 4r^W^ = F a (5.26) 
dG a 

which goes back into Eq.f l5.19l) giving the result 

c a/3 = TY(fr a F (3 e AW = ^2g a5 as/3 . (5.27) 

Putting Eq. fl5.27p into Eq.f l5.23p and using Eq. (15.21 p gives the very simple result 

K^ = -a a0 . (5.28) 
This result satisfies the FDR discussed below in a trivial way: 

Kb = K b p = -W?) (5-29) 

while 

KSl = -^-ImKSl = . (5.30) 



C. Collective part of the vertex function at Second Order 

The collective part K of the two-point vertex function T, at second order is denoted as 
and is obtained by determining fa \ defined as 

0Lr a 

at second order. We have 

Aw{2) = lY, F « F f> G °<p ( 5 - 32 ) 

a/3 

and 

/38 

We then have the standard functional manipulations, 
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5 5G 1 

G a p = — G au r> G V f3 = — G av Gp v Y vrr/ (5.34) 



OLr$ * — ' OCx,? 



where in writing the last equation we introduce the three-point vertex 



= TTT'^aP ■ (5.35) 



5G^ 



Putting this back into Eq. fl5.33p gives 

5AW^ 



^ 2) = -J^— = -ly2 F ^Gu K G, 5 T K5a . (5.36) 

UflKO 

This in turn goes into Eq.f l5.19j) and 



^ _ 

-a/3 2 



G r]g G 7T gT t 



where 



g a p & = Tr (p a (pi3(pse AW (5.37) 

is a three-point self-correlation involving the products of three single particle quantities (J)a\ 
0? and 0^. The second order contributions to the two-point vertex function Y a p can be 
written more symmetrically in terms of the three-point self- vertex 7^5. The latter is defined 

as 



Ga/3-y = — ^ Gaa'Gpp'GsS'la'P'S' ■ (5.38) 
a'P'5' 

Then the collective part of the two-point-vertex is given at second-order by 



K^p — — - JauGuu'G t i t i'G88''lu'^8'C r ^Ti'0'8TT'G v ' v G 7T i 7T T/3 r)7T 
ui/'fi/i'SS' rjrj'mr' 

= —~ ^ lau'S'Guu'GsS'^wq^STvG^'G^iY^^ii 

i/v'/i/i'SS' rjrj'mr' 

o lausG u ^Gs ri Y ' m p (5.39) 



2 
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where in writing the last equality we have introduced the self-dressed propagator 



G a p = } j GqyCr v sGsp ■ (5.40) 

In Appendix O we demonstrate that G a p is a FDRMP with the same generic properties as 
G a p- 

From its definition the collective part of the two-point vertex can be constructed to be 
symmetric and the second-order contribution can be written in the symmetric form: 



Kafi ~ ~ 2 [lausG ur) Gse^rie/3 + JpvsGvrjGsoTrjQal 

— — — llausGu-qGse^-qep + F ai/ sG vv Gso r Yr 1 0p\ 

= — - + r y a us)G uri Gse(T v ep + l^op) — ^ avsG uri G sqT r 



vSr]9 

lausGurjGseJnep , (5-41) 



We look at the properties of K^l in detail in section IVll 



D. Single-particle contribution to vertex function 

Besides the collective contribution to the vertex function r^, we have the single-particle 
contribution. Following the discussion in MMS[7| we find that the single-particle contribu- 
tion 7 a/3 to the memory function T a p at the second order is given by 

= ~ IwsGvnGselnOp ■ (5.42) 

v8r)8 

To avoid ambiguity the collective contribution obtained in eqn. ( I5.4ip is to be denoted 
as K^l. The vertex function is defined by Eq.f l5.38j) and two-point single-particle 
propagator Q is defined in (I5.18p . The quantity G a p is another dressed propagator similar 
to G a/ 3 defined in eqn. (I5.40p corresponding to the collective contribution K a p. G a p is 
defined as 



(XV 

G SrjO tjttQ 'it P 

(5.43) 

vSiyrr 
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In Appendix Owe show that like the set {G^°\G, Q,G}, the matrix G also satisfies the 
property of being a FDRMP. 



VI. ONE-LOOP SELF-ENERGY 



A. General Structure of One-loop Self-energies 

Generally the self energy matrix E is defined in terms of the Dyson equation involving 
inverse of the matrix of the two-point function as follows 

r = G- 1 = G 1 - E . (6.1) 

In the above equation Gq 1 refers to the zeroth order vertex functions. In field theoretic 
models this zeroth order theory stands for the gaussian level theories which usually refers to 
the trivial case e.g., linear dynamics of fluctuations. In the present context the zeroth order ( 
in interaction potentials ) theory refers to the non interacting system. The self energy E will 
be therefore defined here by subtracting the zeroth order and the first order contributions 
from the full vertex function 



- e = r - r<°> - r« . (6.2) 

At the one-loop level the self energy is a sum of collective and single particle contributions 
respectively denoted by K and 7, 

- E = 7 < 2 > + + 0(3) , (6.3) 

where the superscript (2) indicates the second order in in the interaction potential and 0(3) 
stands for higher order contributions. The single particle and the collective contributions 
are respectively expressed as, 



7 (2) 
K (2) 



1 



"I 7 



(G + g)(g + G)-GG-gg 



7 



(6.4) 



-- (jGGT + TGG-f) 
_1( 7 + r)GG( 7 + T) + - a1 GG 1 + ^TGGT 



(6.5) 
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In order to further analyze the one-loop results for 7 and K obtained above we introduce 
the generalized functional Tl[E, A] of a two-point propagator E and three-point vertex A as 

Tl[A,E] = -^EAAE . (6.6) 

The one-loop contributions to the different elements of the self energy E matrix defined in 
eqn. (16. 3p are now expressed in terms of those of matrix II. We focus on the a/i-th element 
(afiE{p, B} ) of the matrix IT: 



Il aiJ ,[A,E}{-l) = - \j ^3^(123)5(1 + 2 + 3)^(2)^(3)^(123) 

K,u,a,S 

(6.7) 

where A is a FDRMP and E is a three-point vertex that satisfies the FDR of section IIII1 In 
the above definition for II, the propagator A is from the set {G^°\ Q, G, G, G} and the three- 
point vertices E include {7^,7, T}. In terms of the functional H[A, E], the single-particle 
and collective contributions respectively given by eqns. (16 .4p and (16.51) are obtained in the 
form: 



7 (2) 



~n [(0 + 0,7] - in [f? )7 ] -in [c, 7 
in [g- 7 + r] - in [g- 7 ] - in [g- r] 



(6.8) 
(6.9) 



Next we consider the fluctuation-dissipation symmetry of the single particle and collective 
contributions to the self enegy S at one loop level. 



B. Partial Self-energies and the FDT 

We investigate the existence of a fluctuation-dissipation theorem (FDT) satisfied by the 
partial self-energy in Fourier space. We first separate n[v4, E] into its components, the 
nonzero choices for the indices k, u, a, and 5. There are nine contributions to n_B_BLA, E], five 
contributions to n p s[A, and five contributions to n_B P [A, E]. In order to better organize 
the algebra we replace all A pp internal lines using the FDR: 
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Ap P (q,u) 



A Bp (q,u) - A pB (q,u) 

i(3u 



(6.10) 



For treating the two-point vertices we introduce a set of simplifying relations which are 
generalizations of eqns. f l3.70p - fl3.72p in terms of the generalized vertex function E : 



"0i 
P 



i(3uiE BiPjPk , 
-i(3 3 u!iU!jUkE 



(3 2 uijUj k E PiBj B k , 



Bi Bj Bk 



(6.11) 
(6.12) 
(6.13) 



where we have explicitly indicated in the first term on LHS the factor of 0, which was earlier 
absorbed in the definition of the frequency u. The quantities fa, ipi and P defined above in 
terms of the generalized vertex function E G {V,^,^ '} satisfies the same set of identities 
given by eqns. (13.731) and (13. 74ft . After considerable algebra we find 

iu x Yl BB [A, E] + U pB [A, E] - U Bp [A, E] 



-— I d2d3 



A Bp (2)A Bp (3)Cl + A pB {2)A pB (3)C, 



- \A Bp (2)A pB (3)C; + A pB (2)A Bp (3)C 2 



5(1 + 2 + 3) 



We have defined the quantities C\, C 2 , and D respectively as, 



(6.14) 



Ci 
C 2 
D 



<Pl{fa + 02 + 03 + 1(>1 + ^2 + 1p3 + P) 
(01 + 02 + ^3) (01 + + 03) - 0203 



(6.15) 
(6.16) 
(6.17) 



Using the identities ( 13. 73ft and ( 13.741) it is straightforward to show that C\ and C 2 both 
vanishes, obtaining a relation similar to (13.341) to hold with the elements of the generalized 
matrix II : 



jr[u] = -ip Ul u B B(i) + n Bp (i) - U pB (l) = 



(6.18) 
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In the above equation and in what follows the functional dependence of T1[E, A] on the 
vertex E and correlation A are not explicitly shown in the RHS to avoid cluttering. In 
Section HTT1 the full two-point vertex funtion matrix was shown to be a FDRMP with the 
result (I3.48p . We have demonstrated here that the self energy matrix TL a p[E, A] expressed 
in terms of the generalized vertex E a ps and correlation functions Aps, are also FDRMP. 
Since T defined in eqn. (I3.35P is linear, using eqn. (I6.18P in the definitions (16. 5p and (16. 4p 
respectively, it follows that 



jr [#(2)] = o , (6.19) 

f [ 7 {2) ] = . (6.20) 

Hence both the collective and single-particle one loop level contributions to the self energy 
respectively denoted by Kij and 7^ are FDRMP. 

After making full use of the vertex identities we find that the self energy matrix element 
ITb p (— 1) is the sum of two pieces 



/ x f d2d3 , 
m Bp {-\) = i / —5(1 + 2 + 3) 



(4>{) 2 A pB (2)A pB (3) + (fo) 2 A pB (2)A Bp (3) 



. (6.21) 



Notice that this self-energy depends only on the three-point vertex functions with one B 
label; E Bpp and E pBp respectively as follows from eqn. (16. lip . 

1. Separation into High and Low Frequency Components 

We will be interested in low-frequency-long-time phenomena. We identify the low- 
frequency contribution to the vertex using the Vertex Theorem (see appendix [B]) as : 

E pBp (0, 0, 0; qi, k 2 , -q% - k 2 ) = Jppp(qi, k 2 , -qi - k 2 ) (6.22) 

where 7 PPP is the static 3-point vertex and write the frequency dependent vertex functions 
as, 



E pBp {ui 1 u 2) uj ? ,]q 1 ,k 2l -qi-k 2 ) = J ppp (qi,k 2 , -qi - k 2 ) 

+ A pBp ( y u 1 ,ui 2 ,ui 3 ;q 1 ,k 2 ,-q 1 - k 2 ) . (6.23) 
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The one loop contribution denoted by T1[E, A] is split in to a low frequency part II (for 
uj— >0) and a high frequency part II 

n Bp (i) = n Bp (i) + n Bp (i) (6.24) 

where n Bp (l) is the contribution with E p b p replaced by its static limit value. After a large 
amount of algebra we find the low frequency contribution 11(1) is reduced to the form, 

rW(i) = ^6A pp {k 2 )A pp {h)[ lppp { qi MM)} 2 (6.25) 

wi A f dx f dy 2 A pp (x,k 2 )A pp (y,k 3 ) 



-oji + x + y — ir] 
where we have defined the operator O as 

~ f d d k 2 f d d k 3 , . . 

The static (oj — > ) contribution 11(1) is given by 

ff Bp (gi, 0) = jp6A pp (k 2 )A pp (k 3 )[ lppp ( qu k 2 , h)} 2 . (6.27) 

Since the low-frequency part II is itself a FDRMP, the BB element of the II matrix is given 
by 



rWi) = im n Bp (i) 

= -^I y ^[7ppp(?i,^2, fc 3 )] 2 A w (x,/c2)^pp(wi -a;, Ala) • (6.28) 
Notice that it follows directly from eqns. (I6.25P and f !6.28|) that the sum rule 

^pfi BB (q,uj) (6.29) 
is satisfied. The high-frequency contribution Hb p (u) to the self energy satisfy 

5^ » Ul (6.30) 
Hb p (1) 
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as Ui — > 0, and hence does not contribute to the slow dynamics of the system. A key 
result in this work is the theorem that Tb pp reduces to the static three-point vertex in the 
low- frequency limit. 



C. Original One loop Problem 

We now focus on the one loop contributions for collective and single-particle self-energies 
respectively denoted by K a/I and 7 aM which were introduced earlier in section PVTl For the col- 
lective contribution K the appropriate three-point vertex functions are E = {(7 + T), 7, T}, 
while the correlation function A = G. The collective contribution is divided in to two parts 
as 



where we have followed the same notations in terms of bar and tilde on the respective terms 
as given in eqn. (I6.24p in indicating the low and high frequency components. Since K ail is 
a FDRMP, we obtain from eqn. (I6.5P and (I6.28p . the low frequency part of the collective 
contribution K^l, as 



For the single particle contribution 7 aA1 the appropriate three-point vertex functions are 



Eqns. (16.321) and (I6.33f) are approximations for the second-order collective and single particle 
contributions to the frequency dependent self energy £ including vertex corrections. We 
write for the respective collective and single particle contributions in the time regime 



K aix {l) = K afl (l) + K afl (l) 



(6.31) 




(6.32) 



E = 7, while the correlation function A = {(Q + G), Q, G}. We obtain from eqn. (16.41) the 
single particle contribution as, 
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Kf B {qxM) = -^OG^t^GppituhhpppTppp • (6-34) 

and 

l^iM) = -^OG pp {t u k 2 )Q pp {t u fc 3 )7 P 2 PP • (6.35) 

In the low- frequency limit we assume, as we show self-consistently, K BB (qi,zi) grows arbi- 
trarily large as the frequency Z\ goes to zero. Since 7^(q , i,0) is regular we can drop the 
single-particle contribution to the memory function in the kinetic equation. 

In the expression given in Eq. (16.341) for the memory function, the implication of having 
the full vertex function T ppp on the dynamic behavior of the fluid has been ignored. To lowest 
order in the interaction potential, the vertex T ppp is approximated by the corresponding 
quantity 7 ppp of the non-interacting theory which is the case for the ideal gas. Since the 
vertex functions are only involved here ( see Appendix B) in the u— >0 limit and hence they 
appears like a static quantity. So far these static approximations are not directly connected 
with HNC or PY approximations. We discuss the role of vertex corrections elsewhere. We 
have then 

r ppp = -7ppp = —jp ■ (6.36) 
We also assume consistent with the long-time approximation that for low-frequencies 



G pp (q,u)^>G pB (q,uj) (6.37) 

and in G pp (q, ui) the G pp (q, u) term dominates and its coefficient can be replaced by its u = 
value. Using the definitions G pp (q,u) = pS(q)F(q,uj) we find 



G pp {q,U)) = (3V{q)pG pp {q,u) 

= V(q)pS(q)F(q,co) , (6.38) 

where we have denoted V = p(3V(q) as the scaled potential. From the above definitions we 

obtain the relation F(q,t) = V(q)S(q)F(q,t). 

We can then write in the long time regime the simple result 
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4 2 i(?, t) = I ^- d F(g - k, t)F(k, t) (6.39) 

We need the Fourier-Laplace transform of Eq. ( 14.331) to go back in the low- frequency form 
of the kinetic equation. After some simple rearrangements we obtain our main result from 
Eq.gSSD 

r^ffe) = [ Me " S Sf p{q - k ' t)p(kA (6 - 40) 

This is a highly nonlinear equation for F(q,t). This result is identical, after reconciling 
notation, to the result in Ref. |6| for Smoluchowsky dynamics. As in SM we can carry out 
an analytic treatment of Eq. (14.331) by expanding 

F(q,z) = ^ + ij(q,z) (6.41) 

for z and lim 2 _>. zip(q, z) — > while \im z ^ ip(q, z) — > oo. f(q) is the nonergodicity 
parameter with the interpretation 

lim G pp (q,t) = f(q)pS(q) . (6.42) 

t— >oc 

In the small z or long time limit the eqn. (16.401) reduces to an intergal equation in terms of 
the nonergodicity parameters f(q) of the form 



1 l! " k k)V(q- k)f(q- k)S{k)V{k)f{k) 



l-f(q) 2p J (2tt) 

(6.43) 

This equation must be supplemented with an equation connecting the potential and the 
static structure factor. At second order in the potential we have 

S-\q) - 1 = V(q) - i J ^~ d V(q - k)S(q - k)V(k)S(k) 

(6.44) 

We choose to fix S(q) and solve for the pseudo-potential V(q). The question of a pseudo- 
potential is discussed in detail in SDENE. We can determine V(q) once we have a form for the 
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static structure factor. Typically we have used the exact solution of the approximate Percus- 



Yevick 



2 if ] equation for hard spheres. One can then carry our the explicit determination of 



the quantities characteristic of the ENE transition. 

The above integral equations for f(q) is solved over a grid to obtain the nonzero solutions 
using iterative methods. The density at which the trivial solution f(q) = changes to a 
nonzero set of nonergodicity parameter values marks the location of the ideal ENE transition 
in the dense liquid approaching from the liquid side. 

In Fig. Owe display the static structure factor for the hard sphere liquid as obtained from 



Percus-Yevick 



211 ] equations with Verlet- Weiss corrections 22j at packing fraction rj = 0.62. 



Using this S(k) as an input we solve iteratively eqns. (16.44j) to obtain the renormalized 
potential V(q). The wave vector q is chosen over a grid of upper cutoff value qa = 80 and 
having 500 points. Fig. O displays the V(q) obtained with input structure factor of Fig. O 
Next we use this renormalized potential to evaluate the vertex functions in eqn. ( I6.43|) for 
the the nonergodicity parameters f(q). Solving these integral equations iteratively we obtain 
that the f(q) are vanishing till the critical packing fraction of r\ = .62. The corresponding 
nonergodicity parameters f(q) over the whole wave vector grid is shown in Fig. O The ENE 
transition point of the core problem of self-consistent feed back mechanism using Percus- 
Yevick strucutre factors with Verlet- Weiss corrections therefore at packing fraction of .62 
which is close to the close-packing density. 



VII. CONCLUSIONS 



We have shown within the pseudo-potential expansion the ENE problem for ND reduces 
to precisely the same problem for SD. For hard-spheres and using the PY approximation 
for the structure factor we can solve the approximate ENE problem to find 77* = 0.62, and 
the two-step exponents a and b consistent with the results found in SM. This is a nontrivial 
result and is closely associated with the vertex theorem which relates the low-frequency limit 
of r_B PP to the static vertex 7 ppp . While we have proved this result only to second order in 
perturbation theory, we expect it is more general [13 1. 

We have found simplification of the ND case only in the low-frquency regime. More 
generally the SD case is considerably simpler as demonstrated in MMS where the zeroth 
order vertices are very simple. 
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This result for the dynamic structure factor raises questions about higher-order correla- 
tion functions such as G ppp (123). There is good reason to believe that we can make progress 
in determining this quantity in the low-frequency regime. 

Our results here depend on ignoring static-vertex corrections and higher-order loop con- 
tributions. We need to investigate these correction to see if the results for the ENE are 
stable. We see from Eq. (16.341) that the low-frequency dynamics depends not just on the 
static structure factor, but also on the static three-point vertex. 

The theory presented in this work has been set up to deal with what we call the core 
problem in the ND case. The primary motivation here is the determination of the observables 
involving these core variables $o — {PiB). The terms core refers to the fact that the 
Hamiltonian and hence the corresponding MSR action is expressed in terms of only these 
variables. As already pointed out ( see sec. 9 in ND1), by going beyond the core variables, 
considerations of additional degrees of freedom enter the theory. For example in the ND 
case, we have a larger phase-space due to the momentum degrees of freedom and there 
are additional conservation laws for the system. If one extends the set of core variables to 
include the momentum density, we have all the correlation functions among p, B and g in the 
description. The resulting hierarchical structure in the formal expressions for the correlation 
functions may change the viability of the ENE transition being discussed here. In this regard 
it is useful to note the known results of including the momentum density fluctuations in the 
fluctuating hydrodynamic approach to the problem. The momentum density coupling to 
the density fluctuations in the hydrodynamic equations give rise to the ergodicity restoring 
mechanisms (23I 24\. In the case of SD on the other hand, the core problem itself covers 
essentially all of the degrees of freedom of interest. However the microscopic dynamics in this 



25 



26) 



case is dissipative and the equations of motion for the collective collective variables 
involve multiplicative noise which may have implications in restoring ergodic behavior in 
a liquid. The present work demonstrates that in the case of the core problem through a 
reorganization of the perturbation theory in terms of an effective potential, the ideal ENE 
transition is pushed to the close pack density for a hard sphere system. 
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Appendix A: Three point vertex functions 



We obtain the FDT relations between the different three-point vertex functions here. For 
this we use the general definition of the three point correlation functions and the correspond- 
ing three-point vertex function as, 

G aiHUk = -G aa >(i)G^r(j)G w r(k)r a y. v ' k (Al) 

where a.i stands for the field a at the point i and so on where everything is in Fourier space. 
The labels a,p and v are respectively taken from the set {p, B}. Using the result that 
^piPjpk = 0' we obtain from the formula flAlj) . the various cumulants are : 

-G PiP . Pk = G pp (i)G pp (j)G pB (k)F PiP . Bk + G pp (i)G pB {j)G pp (k)F PiBjPk 

+ G pB {i)G pp {j)G pp (k) + G pp {i)G pB (j)G pB (k)T PiBjBk 

+ G pB (i)G pp (j)G pB (k)Y BiPjBk + G pB (i)G pB (j)G pp (k) 

BiBjp k 

+ G pB (i)G pB {j)G pB {k)T BiBjBk (A2) 

-G BlPjPk = G Bp (i)G pB {j)G pp {k)Y p . B . Pk + G Bp (i)G pp (j)G pB (k)T PiPjBk 

+ G Bp (i)G pB (j)G pB (k)T PiBjBk (A3) 

—G BiBjPk = G Bp (i)G Bp (j)G pB (k)r PiP . Bk . (A4) 

Substituting the results flA2j) - flA4j) in eqn. (I3.54p . we obtain an expansion in terms of the 
various three-point vertex functions. For example the coefficient of the vertex function 
^B iP]Pk is obtained after some trivial but tedious algebra as 

^B iPjPk ■ mG pB (i)G Bp (j)G Bp (k) +cuiU)jG pB (i)G pp (j)G Bp (k) 

= i^iG pB {i)G pB {j)^G Bp {k) - iu k G pp (k)j = iutG pB {i)G pB (j)G pB {k) (A5) 

In reaching the above result we have used the FDT relation (I3.48p . Similarly for the other 
two vertices with one B field, i.e., T p . B . Pk and T PiP:jBk are obtained. For the vertex with two 
B fields, e.g, T BiBjPk we obtain the corresponding coefficient as 



^BiB jPk ■ UiUjG pB (i)G pB (j)G Bp (k) - iuiUjUkG pB (i)G pB {j)G pp (k) 
= UiUjG pB (i)G pB {j)G pB (k) 
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(A6) 



and similarly we obtain the coefficients of the other two vertices each with two B fields. 
After organizing the coefficients of the different three-point vertex functions the following 
result is obtained. 

i{uiT BiPjPk + u}jT PiBjPk + u) k r p .p.Bk) - iUiUJjOJ k r BiBjBk 
+ \u i u; j r Bi B j p k + UjUj k T PiBjBk + u k UiT BiPjBk } = (A7) 

In reaching the above result we have dropped a nonzero common factor of 
G pB {i)G pB {j)G pB {k) from the LHS. 

Next, consider the relations (I3.55P involving three-point cumulants having two B fields. 
Once again substituting the results (IA3p -f lA"4|) in eqn. (I3.55p . we obtain after organizing the 
coefficients of the different three-point vertex functions the following result. 

G Bp {k)^iu k Y* p . p . Bk +iujiY BiPjPh + iujjT PiBjPk +UiUjT Bi B jPk j 
= G pB (k)^iuiT BiPjPk + iu}jT PiBjPk + iu k T PtP]Bk 

+ UiU)jr BiBjPk + UiU k T BiPjBk + UjU k T PiBjBk - toJiWjUkTiuBjBi } = . (A8) 

In obtaining the above result, we have dropped the common factor of G pB (i)G pB (j) from 
both sides and used the basic FD relation (I3.48p . Using eqn. f!3.65|) in the result (IA8I) we 
find that the coefficient of G pB {k) vanishes and one has the result 

iuJ kT* PiP]Bk + iuJiT BtPjPk + /-,r ( , /;,./, + 1 '/.•/;,/,. = . (A9) 

Appendix B: Higher-Order Thermodynamic Sum Rule 
1. Three-point Quantities 

Let us consider the low frequency behavior of the full three-point cumulant G pBB (123). 
The quantity p(l) in G pBB (123) can be replaced by an arbitrary function of density as long 
as each density corresponds to the same time. The three-point vertex r Spp (123) is related 
to the 3-point correlation G pBB by the general relation 
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r Bpp (i23) = -r Bp (i)r pB (2)r>(3)G pBB (i23) . (bi) 

Among the FDR identities for the 3-point cumulants, we have for the imaginary part of 
GpBB the relation : 

G^(123) = ^G; pB (123) + ^G; Bp (123) . (B2) 
Next, look at the definition of the inverse time fourier transform 

G^bWMMM) =f^f ^e-^^-^e-^^-^G^BCwi,^) (B3) 
J 2ir J In 

where we have introduced 

G Q/ 3 7 (123) = G a/3l (u)i,uj2, qi, q2, qz)S(ui + u 2 + u 3 ) . (B4) 
In this notation we do not always write the third frequency entry since it is implied: 

G a/ 3~,(u)i, uj 2 ) = G a fa(ui, lu 2 , -oji - u) 2 ) (B5) 



and the wave number dependence has been suppressed. Setting t 2 = t 3 in eqn. (IB 3 1) we 
obtain, 

G pBB (t 1 ,t 2 ,t 2 ) = J^J ^e-^-^GpBBiuu^) . (B6) 
This vanishes for t 2 > t\. This is consistent with 

^■G p bb(ui,u} 2 ) (B7) 

being analytic in the UHP for oj\. Assuming G p bb(^i,^2) is analytic in the UHP for oj± we 
can write a dispersion relation 

C,»K«*) = f^S^L . (B8) 
Putting Eq.flBU) in Eq.(jB2]) gives 

G" pB B^2) = ^GUk,^) - ^f^ G'^uuUi) . (B9) 
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Putting Eq.flBlD in Eq.flBg) gives 

d(D /3u 2 G' ppB (q; Q, u 2 ) - p[Q + u 2 }G' pBp (u, u 2 ) 



G p bb{^Xi u 2 ) — J 



2tt u) — u>i — irj 

Letting lo\ and u 2 go to zero gives 

G pB b(0,0) = -/3 J ^G' pBp (u,0) 

We have a FDR identity 



which tells us that 



and we can write 



£g pBp (u,0) . 



In the time domain 



G pBp {tiMM = / ^e-- (tl - t3) e--^-^G pBp (o; 1 ,c 2 ) 

Fourier transforming over t 2 obtains 

G pBp (h,u 2 ,t 3 )= / ^e-^-^e^GpB^u^-ut - u 2 ) 



2tt 

Setting t 3 = ti and letting uj 2 — >■ 0: 



"2 — GpBp{tOi : 0, — wi) 

Combining Eq.f ]B14| ) and ( 1B17|) gives 



G p bb(0,0,0) = -/3G p bp(*i,0,*i) 
When the times of the p's are equal in G pBp we have 
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G P Bp(q;ti,t 2 ,t 1 ) = J^I ^e-^-^G^^u,) (B19) 
Introducing f(ti) = Sp(qiti)Sp(q 2 ,ti) one has Gf B (0) satisfies the two-time FDR 

G fB (0) = -fij ^G fp {u) = -(3(5p(q 3 )f) (B20) 
= -PS ppp ( qi ,q 2 ,q 3 ) (B21) 
and the static three-point cumulant enters the development. We have then 

G pBB (q; 0, 0) = -(3G fB (q; 0) = f3 2 S ppp ( qi , q 2 , q 3 ) . (B22) 
In terms of the three-point vertex 



r Bpp (0,0,0) = -T Bp (0)T pB (0)T pB (0)G pBB (0,0,0) 

P 2 S ppp (qi,q 2 ,q 3 ) 



(-PS pp ( qi ))(-pS pp (q 2 ))(-pS pp (q 3 )) 
= -P~\ PP (qi,q2,qs) • (B23) 

which is a result of much use in evaluating the one-loop contribution to the self energy both 
at single-particle and collective levels. Note that 7 PPP is a static three-point vertex. 

Appendix C: FDR Matrix Propagators 

FDR matrix propagators (FDRMP) A pi/ (q,u) satisfy the following properties: 

A^(q,uj) = A* (q,u) 



+i(3uA pp (q,oj) = A Bp (q,u) - A pB (q,u) 

dCu (3cuA pp (q, Co) 



A pB (q,u) = J 



2n ui — u + irj 



(CI) 



From this it follows 



A pB (q,0) = - J ^-(3A pp (q,u) . (C2) 
Finally the element A BB (q,u) = 0. 
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We now prove the following important property of the FDRMP : If A a/ 3(q,u) and 
C a p(q,uj) are FDR matrix propagators then 

D a/3 (q, u) = ^ A <*n{q, ^)^^u(q)C^(q, u) (C3) 

is also a FDR matrix propagator. The proof is rather direct. Look first at the response 
channel: 



D BB (q,u) = ^2A Bf ,(q,u)a Bu (q)C u(3 (q,uj) 

= A Bp (q,uj)a pp (q)C pB (q,uj) = . (C4) 
Consider next the off -diagonal components 



D pB (q,u) = A pB (q,u)a Bp (q)C pB (q,u) = A pB (q,u)V{q)C pB (q,u) 

D Bp (q,u) = A Bp (q,u)V(q)C Bp (q,Uj) (C5) 

It is easy to see that D pB (q,u) = D B (q,u>). Next consider the diagonal component 



D pp (q,u) = A pp (q,u)V(q)C Bp (q,u) + A pB (q,u)V(q)C pp (q,u 
V(q) 



i/3uj 
V(q) 



[(A Bp (q,u) - A pB (q,u))C Bp {q,Uj) + A pB (q,uj)(C Bp (q,Uj) - C pB (q,u))\ 
[(A Bp (q, u)C Bp (q, u) - A pB (q, u)C pB (q, w))] 



i(3u 

The above result implies that 



' [D Bp (q,u)-D pB (q,uj)] (C6) 



i/3u)D pp (q,cj) = D Bp (q, u) - D pB (q, uj) (C7) 
-tf3L0D* pp (q, u) = D pB (q, u) - D Bp (q 1 u) (C8) 

Together Eqs.tfUTD and (EI]) give D pp (q,u) = D* p (q,u). 

Let us now consider the dressed propagators respectively denoted as G and G. In operator 
notation G and G are respectively defined as G = QoG and G = QaGaQ. Writing out 
explicitly the matrix forms we obtain for G and G the following expressions : 
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G a/3 (q,uj) = G afl (q,uj)(r fMU (q)G vl 3(q,uj) 

G a/3 (q,uj) = G a ^(q,Uj)(T^ u (q)G uS (q,Uj)a Sv (q)g vl3 (q,uj) 



(C9) 
(CIO) 



From the above theorem then it follows that both G and G are FDRMP. This holds since 
each of the G, Q, and G^ satisfies the conditions of being a FDRMP. 
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FIG. 2: The renormalized potential V(q) corresponding to the structure factor displayed in Fig. 
1. The inset shows the oscillations at large wave vectors on a enlarged scale. 
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qo 



FIG. 3: The non ergodicity parameter f(q) vs. qo at packing fraction r] = .62. The upper inset 
shows the same quantity at r/ = .60 indicating ergodic behavior. The lower inset shows the f(q) 
vs. qo at the ENE transition point rj = .525 in the earlier MCT of Ref. 10 
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